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I. INTRODUCTION 



In recent years, theoretical and computational studies of molecular spectra have been 
one of the most valuable tools available for studying on atoms and molecules. At its sim- 
plest level, knowledge of spectral characteristics allows us to detect the presence of particular 
characteristic and essential components of matter. Especially, molecular spectra can be used 
to understand the motion of electrons in molecules as well as the vibration and rotation of 
the nuclei. The chemical interactions between atoms and molecules assist to investigate 
the physical properties of individual molecules. In the light of this knowledge, dissociation 
channels 111, centrifugal distortion constants t2|, semiempirical dipole moment functions tSl 

nnriQ 

and other data about the rotation, vibration and electronic energy levels 1^, |5|, |6|, [TJ of di- 
atomic molecules can be accurately determined by using theoretical methods. Moreover, 
some quantum-mechanical calculations on rotational and vibrational energy levels of di- 
atomic molecules have been applied to problems in molecular physics for a number of years 
Sj. The modified shifted large 1/N approach has been applied to obtain energy levels of a 
rotational potential 9|, arbitrary £-state solutions of the rotating Morse potential has been 
investigated through the exact quantization rule method [lO| and other algebraic approaches 
and applications have 



potentials 
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jeen previously applied to rotational and vibrational states of rotating 



In this study, the bound state energ y le vels are obtained by solving the Schrodinger 
equation for the Morse jl8] and Kratzer jiol l molecular potentials together with an exactly 



solvable angle-dependent potential, respectively. 



Veie) 



(1) 



r — r.- 



Ve{9) 



(2) 



where subscripts M and K indicate the Morse and Kratzer potentials, respectively, a con- 
trols the width of the potential and is the equilibrium internuclear distance. The quantity 
De is the electronic (or spectroscopic) dissociation energy of the diatomic molecule and it 
differs slightly from the chemical dissociation energy Dq, i.e., Dq = De — huJe/2, where Ue is 
called harmonic vibrational parameter 20, l2l|]. Moreover, the minimum value of VM,K{r) at 
r = belongs to Df.. The second term in the right-hand side of Eq.([T]) or Eq.([2]) represents 
an angle-dependent potential and its uncovered form is given as follows 



A and B in Eq.Q are fixed constants or parameters obtained by some fitting procedure 
which is based on experimental or theoretical results; it is important to emphasize that they 
cannot depend on the angle 9. The factor /i^/2/i is introduced in view of future convenience. 
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The potential given in Eq.Q has been introduced for the first time by Makarov et.al 
classifying some non-central potential systems. The Schrodinger equation for this type of 
angle-dependent potential can be exactly solved to obtain the bound state energies of a 
diatomic molecule. It is well-known that the problem of exact solution of the Schrodinger 
equation for a number of special potentials has been a line of great interest in some quantum 
mechanical applications. The solution of this equation for some potential has been made 
by applying some analytical methods. One of these methods is developed by Nikiforov 
and Uvarov [23] as a new approach to the theory of special functions. They succeeded 
in obtaining an unified integral representation for functions of hypergeometric type. This 
type of hypergeometric equation with an appropriate coordinate transformation is given as 
follows 



r{s) + ^^'(s) + ^^(.) = (4) 

where a{s) and a (s) are polynomials, at most second— degree, and r (s) is a 
first— degree polynomial. The general view point of this paper is to present an analytical so- 
lution of the angle-dependent part of the Schrodinger equation for an exactly solvable angle- 
dependent potential Ve(6')/r^ and also to obtain modified i states of diatomic molecules. 
The solution method developed by Nikiforov and Uvarov is used for solving the Schrodinger 
equation. The angle-dependent part of the Schrodinger equation is investigated in detail 
to derive some analytical result and the solution of the radial part of the associated equa- 
;ion for the Morse and Kratzer potentials is extracted from the papers published previously 
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251]. The modified i state expressions for the Morse and Kratzer potentials are obtained 
by connecting the results of the angle-dependent part with the radial one. The modified i 
states of a number of neutral diatomic molecules composed of a first-row transition metal 
and main-group elements are calculated for both Morse and Kratzer potentials with an 
angle-dependent potential. 
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II. SEPARATING VARIABLES OF THE SCHRODINGER EQUATION IN 
SPHERICAL COORDINATES 



The starting point of this section is to separate the Schrodinger equation in spherical 
coordinates for a diatomic molecule represented by a rotating potential model. After sepa- 
rating the center of mass motion, the eigenvalue equation for a rotating motion in spherical 
coordinates is solved by using the NU method and the energy levels of the discrete spec- 
trum are obtained for several diatomic molecules. In spherical coordinates, the Schrodinger 
equation is written as follows: 



2^ 



1 d 



y.2 Qy, \ Qy, 



d 



+ ^ ^ (sine— 
r'^sinO 39 \ 39 



1 



3^ 



9 ■ 9/1 o 9 I ^ nlmiy) 

r'^sin 9 3(p j 



(5) 



The energy E in Eq.([5]) is real and it is either discrete for bound states {E < 0) or 
continuous for scattering states {E > 0). Introducing a new variable x = cos'^9, Eq.Q can 
be explicitly turned into the more useful one: 



1 3 
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J.2 Qf \ Qr 



3'^ 3 13'^ 
4x(l - + 2(1 - 3x)^ + 

3x^ 3x 1 — X a(y9^ 

2/i 



^n£m(r) 

0. 



+ ^{E - V(r))^„,„(r) = 
Consequently, this equation is separable for a potential of the following form. 



(6) 



^(r) = VM,K{r) + \ 



Veix) + 



1 



1 — X 



(7) 



If we write the wave function as \E'„£m(r) = r^^Rne{r)Q£m{9)^m{'^), then the wave equa- 
tion in Eq.([6]) with the potential in Eq.Q is separated to a set of second-order differential 
equations in all three coordinates as follows: 



Eg ^ 2l2 



V) = 0, 



4x(l - + 2(1 - 3x)^ - ^ + - ^^Ve{x) ) e,^(x) = 0, 



dip"^ 



'^V^{^) + E^ $™(^) = 0, 



(8) 

(9) 

(10) 
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where E^p and Eg are the separation constants, which are real and dimensionless. Since 
the wave function \l/„^m(r) must be finite in all space for the bound states, the boundary 
conditions for Eq.([H]) require RntiS^) = and the square-integrability of Rniij) on (0, oo), 
which implies that i?„^(oo) = 0. The finite solutions for Q^mid) in the range < < vr are 
able to map into a differential equation of hypergeometric type. Moreover, the boundary 
conditions for Eq. (fTO|) must be ^m{,V + 2vr) = If the azimuthal-dependent potential 

part Vip{ip) is set up to zero, then the normalized solution of Eq. ffTOj) that satisfies the 
boundary conditions becomes 

$^(¥;) = ^e™^ m = 0,±l,±2,..., (11) 
V Svr 

where one of the separation constants E^^ represents m^, i.e., E^^ = m^. 



A. The Solution of Eq.® 

It is well-known that the solution of the radial part of the Schrodinger equation gives 
eigenvalues and eigenf unctions for a particle moving within the interaction potentials. How- 
ever, the solution of the angle- dependent part of the corresponding equation does not de- 
pends on eigenvalues presented in the solution of the radial part explicitly. It only exhibits 
a parameter relationship between contribution constants which come from the ^-dependent 
part of the potential. Such a relationship can be expressed by solving Eq.([n]) in terms of 
Eg. Eq.Q can then be rewritten in the following form by introducing an exactly solvable 
angle-dependent potential given in Eq.Q, 



(ax{1 - x)^ + 2(1 - 3x)-^ - + Ee~ ( + -]] QU^) = 0. (12) 

\ ax"^ ax 1 — X \1 — a; x J J 

An arrangement of the above equation turns to a convenient form to make a comparison 
with the main equation of the NU method given in Eq.(j4]); 



(P&em{x) (1 - 3x) d&im{x) 



dx"^ 



2x(l - x) dx [2x(l - x)f 
(-Eex" + x{Ee -A + B)-B^ Qi^{x) = 0, 



(13) 



5 



where A = m? + A (keeping in mind the selection of i?,^ = m^). Having compared Eq. f|T3|) 
with Eq.(jlj), the following polynomial equalities are obtained immediately 

r = l-3x, (14) 
(T = 2x(l-x), (15) 



cr 



-Eex^ + x{Ee -A + B)-B. (16) 

In the next step, the basic solution procedure of the NU method given in Ref. will be 
followed to find a solution of Ea. ffT3l) in terms of Eg. If polynomials given in Eqs.(|T4l)-( |T6ll 
are substituted into Eq.(6) of Ref. 2J], tt function is obtained as follows 



TT = ± -\jx\AEe - 8fc + 1) - x{AEe - 4A + 45 - Sfc + 2) + 1 + 45. (17) 

The simplest form of tt can be written 

vr = ^^±^Vax2-/3a; + 7, (18) 

where a = AEg -8k + l, (3 = AEg - AA + AB-8k + 2 and 7 = 1 + 45. The possible solutions 
according to the plus and minus signs of Eq. f|T8l) depend on the parameter k within the square 
root sign. The expression under the square root has to be the square of a polynomial, since 
TT is a polynomial of degree at most 1. To satisfy this condition, the discriminant of the 
expression within the square root must be set up to zero, i.e., A = — 4a7 = 0. This 
identity leads to 

{AEe-AA + AB-8k + 2Y~A{AEg-8k + l){l + AB) = 0, (19) 
and a second-order equation related to k is originated as follows 

4fc2 + Ak{A + B-Eg) + {A- Bf - 2Eg{A + 5) + E^ 
Hence, the double roots of k are derived as 



- A = 0. (20) 
(21) 
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Substituting fci 2 into Eq. (fT7|) . the four possible solutions of vr are obtained 



TC 



— X 1 



for = -^^±^ + y A{1 + 4B) 



(22) 



2VA + VI + 45 j X - VI + 45j , 
for fc, = - (^+^^-^^) -iJl(l + 45) 



where ki^2 is determined by means of the same procedure as in Ref. 2J|. We have to 
choose one of the four possible forms of tt to obtain the bound state solutions. Therefore, 

for 



its most suitable form is established by vr = 



2VA + VI + 45 X - VI +45 



k 



(A+B-Eg) 



2 2 v + 45). The main requirement in the selection of this form is to 

find the negative derivative of r(s) given by Eq.(9) of Ref. [2^. In that case, r(s) and t'{s) 



are obtained, respectively. 



TiS 



1 + Vl + 45 -x(4 + + VTT 



45 



r'is] 



- (4 + 2VA + Vl + 45) < . 



(23) 



Another major polynomials given in the basic solution procedure of the NU method are 
A and An 23|. Both polynomials can be connected with each other by means of Eq.(7) and 
Eq.(8) of Ref. 24 1. Hence, a polynomial of degree n is found by using \n = —nr' — "^""'"'V ^^; 



\n = + 2n + 2nVA + n Vl + 45, (n = 0, 1, 2, ...) 
taking a" = —4. Moreover, A is obtained from k2 + tt'; 

A = _ivrT45 (l + -^(l + 5-5e + l) -Va 



(24) 



(25) 



After comparing Eq.( ]24l) with Eq.( ]25ll and also making some arrangements on the com- 
parison, the separation constant Eg is obtained as follows 



2n+VA) + 2VA + Vl + 45 + 2n + V A) Vl + 45 + (1 + 5) = Eg. (26) 



It is very useful to prepare i{i + 1) as a new presentation instead of Eg. In this case, 
Eq.(|26l) turns to 



(l/2 + 2n+^+ Vl/4 + {l/2 + 2n+^+ ^1/4 + 5 + l) =£(£ + 1), (27) 



and it becomes in terms of 



(^1/2 + 2n + + V^TTT 



B 



(28) 



The term I in Eg. (128!) can be named the "modified" orbital angular momentum, since 
the contribution which comes from the angle- dependent potential damages the usual orbital 
angular momentum £. Moreover, the result obtained in Eq.( |28l) is in agreement with results 
on the more involved case of Ref. 26|. In the limiting case 5 = 0, the factor a/1/4 + 5 in 



Eg. (1281) should be replaced by ±1/2 so that Eq.( |28l) turns into v+ V A, where v = l + 2n for 



271]. The parameter 



the odd functional solution oi v = 2n for the even functional solution 
I does not need to be integer. However, the difference between the parameter £ and the 
square root terms in Eq. (!28|) have to be integer; 



(29) 



n = i{^- (l/2± V^± v/l/4±5)}, n = 0,l,2, 
where n corresponds to the number of quanta for oscillations. 



B. The Solution of Eq.® 

It is remarkable that the radial equation in Eq. ([H]) is independent of the angle-dependent 
term given in Eqs.([T]) and ([2]) for the Morse and Kratzer cases, respectively. Eq.([8]) is exactly 
soluble by means of the NU method. However, some caution must be observed especially 
on the solution of the Morse potential since the exponential nature of the Morse potential 
and the radial behavior of the centrifugal kinetic energy term do not allow for solving 
the Schrodinger equation simultaneously. In the case of Kratzer potential, no caution is 
necessary when considering the Kratzer potential together with the centrifugal term since 
both terms shows the radial behaviors. In the following subsections, the solution of both 
potentials is briefly investigated by using the NU method. 
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1. The Morse case 



Adopting the Morse potential to Eq.®, the radial Schrodinger equation turns into the 
following form 



+ '^{E-D, {e-^<--'^) - 2e-"('--^'=)))^ i?„,(r) = 0. (30) 

Disadvantage of Eq. ( !30l) is that analytical solutions cannot be found because of the centrifu- 
gal kinetic energy term of the potential proportional to Ee/r"^ is included into the radial 
Schrodinger equation. In order to obtain an analytical solution of Eq.( l30l) . the term Eg/r'^ 
has to be approximated to the exponential one. Using an accurate approximate treatment 
suggested by Pekeris 28|, this term can be translated into the following form 



Ee Eg , 



where x is a coordinate transformation represented by {r — rf.)/rf. and Di is the coefficients 
which are given in Eq.(18) of Ref . |24| (i=0,l,2). Substituting Eq. (l3T]) into Eq. fl30|) and using 



a new variable of the form s 



-ar^x 



the resulting Schrodinger equation becomes 



d^Rni{s) 



1 dRrifis] 1 r n 01 -r, / N 



(32) 



where -el = {e - /h'a^ = 2/i (2De - ^) and £3 = 

2/i (^De — -^fl^ j /h'^a^. Comparing this equation with that of Eq.(21) of Ref. 24| and follow- 
ing the solution steps of the NU method, the energy spectrum according to the quantum 
numbers ra, n and m is obtained as 



where 



E ~ 



n^Eg 

2/xr2 



1 - 



arp 



2^2 



h a 



1 

n -\ — 
2 



2^ 



Eg_ /J_ 

r1 KaTf. 



(33) 



(34) 



and Eg is given by Eq.( l27|) . keeping in mind A = m? + A. The highest vibrational quantum 
number Umax can be directly estimated from the condition dEnnm/dn = 0; 

nmax is generally limited to obtain the number of bound states in the case of the Morse 
potential and its maximum value depends on the potential parameters of a given diatomic 
molecule as well as the quantum numbers n and m. 



2. The Kratzer case 



Among many two-particle interaction models, one of the most interesting potential types 
is the Kratzer potential because it can be exactly solved for the general case of rotation 
states different from zero. The first term on the right-hand side of Eq.Q is the central 
Kratzer potential and the radial part of the Schrodinger equation in the presence of this 
potential can be written as follows, recalling Eq.Q, 

2" 



2i2 



E + De-De 



0. 



(36) 



Using the transformation s ^ r/r^ and letting the dimensionless notations 



-el 



2firlE 



2 ' 



e2 



e^ = Ee 



(37) 



Eq. (l36!l can be rewritten in a simple form as follows 



d Rnti^s) 1/ 22, \ TD / \ 

^^2 + ^ (-^1^ + - Rni{s) 



0. 



(38) 



The complete solution of Eq. fl38l) by means of the NU method can be found in Ref. 25| . 
after having made of some notation setting. Hence, the energy spectrum with respect to the 
quantum numbers n, n and m is obtained as 



Ennm 



'212 



l + 2n+ Vl + ADr, 



(39) 



where 



Dr. 



X 



+ (1/2 + 2n + Vm^ + A + ^1/4 + 5) 
(1/2 + 2n + Vm^ + A + ^1/4 + B + 1 



(40) 



The derivative of Eq.( l40l) according to n gives the maximum vibrational quantum number 
Umax in the case of Kratzer potential; 



dE ~ 



0. 



(41) 



dn {l + 2n^a. + VT+W^y 

The condition which requires to satisfy the equality on the right-hand side of Eq. (l4Tl) is that 
Umax niust bc supportcd by an infinite number of vibrational levels. 
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C. Remarks and calculations for the modified ^ states 



In order to discuss the behavior of energy spectrums of a diatomic molecule when the 
values of quantum numbers n, n and m differ, it is very useful to select some diatomic 
molecules composed of a first-row transition metal and main-group elements (H-F). One or 
two of these molecules are the first-row transition metal hydrides such as ScH, TiH, VH, 
CrH and MnH |29]. Transition metal hydrides are chemical compounds formed when hy- 
drogen gas reacts with transition metal atoms. These are of considerable importance in 
chemical synthesis as intermediates and in solid matrix samples for infrared spectroscopic 
study. Another diatomic molecule containing the transition metal element copper (Cu) and 
the main group element lithium (Li) is CuLi, which elucidates the nature of the bonding in 



mixed transition metal lithides 



301 1 . Presently the transition metal carbide molecules such 



as TiC and NiC represent a very active field of research, especially due to the desire for a 



quantitative understanding of their chemical bonds 



31 



32j . Moreover, diatomic scandium 



nitride molecule ScN has excellent physical properties of high temperature stability as well 



as electronic transport properties, which are typical of transition metal nitride 



Fur- 



thermore, the scandium fluoride molecule ScF is the best studied transition metal hahde 
and it has been fairly well characterized 3J| . Diatomic molecules which consist of transition 
metal and main group elements are challenging theoretically and computationally, but re- 
cent advancements in computational methods have made such molecules more accessible to 
investigations. Their spectroscopic parameters have been accurately determined by using ab- 
initio calculations. One of these calculations is called the multi-configuration self-consistent 
field (MCSCF) and it seems qualitatively correct. In Table 1, the spectroscopic parameters 



of the above mentioned diatomic molecules are summarized using MCSCF results |35|. How- 
ever, choice of the parameter a is not a simple issue. Solution of the Schrodinger equation 
for the Morse potential gives the following well-known relation (see p. 132 of Ref. 2l|): 

(42) 



where B^. = g^2^e^^2 • Notice that the parameter a is used to calculate the energy spectrum of 
the Morse potential. Another considerable effort for the Morse potential is that the highest 
vibrational quantum number Umax changes according to the spectroscopic parameters of 
diatomic molecules as well as the parameters n, m, A and B, keeping in mind Eq. (l35l) . As 
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an example, the value of Umax for ScH is 20 in the fixed values of A = 1 and B = 9 and 
under the conditions of < 10 and m < 10. The values of rimax for TiH, VH, CrH, MnH, 
CuLi, TiC, NiC, ScN and ScF molecules given in Table 1 are aligned 20, 20, 17, 14, 70, 71, 
50, 100 and 131, respectively, in the same values of parameters and conditions. 

To calculate the bound state energies of diatomic molecules given in Table 1, Eqs. fl33p 
and (139!) must be recalled for the Morse and Kratzer cases, respectively. Taken into account 
spectroscopic parameters of diatomic molecules and arbitrary values of A and B, the bound 
state energies can be compared for both potentials. This type of comparison is given in 
Table 2. As can be seen from Table 2, when parameters A and B are fixed to 1 and 1, 
respectively, for different values of n, n and m, the bound state energies become lower than 
that of other values. A comparison of A = 1 and B = 9 with A = 9 and B = 1 shows that 
the bound state energies obtained for A = 1 and B = 9 are a little smaller than the energies 
obtained for A = 9 and i? = 1 in small values of n, n and m, especially and 1. For large 
values of the quantum numbers, the bound state energies obtained for A = 9 and B = 1 
tend to become more separately spaced than the energies obtained for A = 1 and B = 9. 



III. CONCLUSIONS 



An interesting extension of this work is to study the effect of an angle-dependent poten- 
tial to the Morse and Kratzer potentials and to examine the partial changes on the usual 
£ states. The analysis presented in this work suggests that the bound state energies of 
diatomic molecules depend on the quantum numbers n, n, m and also the parameters A 
and B. Moreover, the energy spectrum obtained in Eq. (!33|) is an approximate description 
of the quantum aspects of diatomic molecules for the Morse potential together with angle- 
dependent potential while the spectrum obtained in Eq. (l39l) is a complete description for the 
Kratzer potential together with angle- dependent potential. Furthermore, the solution proce- 
dure presented in this paper is also systematical and efficient for solving the angle-dependent 
part of the Schrodinger equation. 
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TABLE I: Spectroscopic parameters and reduced masses for some diatomic molecules composed of 
a first-row transition metal and main-group elements (H-F). The complete list of this table can be 
found from Ref. 35 1. 



Molecule (eV) (A) LOe {cm ^) a (A ^) fi (a.m.u) Reference 
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\m 


MnH 
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m 


CuLi 


1.74 


2.310 
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1.00818 
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m 


TiC 


2.66 


1.790 
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9.606079 


m 


NiC 


2.76 


1.621 
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9.974265 


[32] 


ScN 


4.56 
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[33] 


ScF 
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16 



varif 
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ttion of bound state energies (in eV) for various values of n, n, m, A and B 



Morse Potential 



Kratzer Potential 



A = l 
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A=l 
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-1.60579 
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-1.56867 
-1.56497 
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-2.54773 
-2.40344 
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-5.72965 
-5.72723 
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a(r-re) 

FIG. 1: Representative vibrational energy levels and rotation of a diatomic molecule, n is the 
vibration quantum number and Dq is the chemical dissociation energy of the lowest (n = 0) 
vibrational level. The internuclear distance r is representatively shown in the right-hand sight of 
figure. 
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